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Abstract 

The applied magnetic field and temperature phase diagram of the archetypal frustrated mag¬ 
net, Gd 3 Ga 50 i 2 , has been reinvestigated using single crystal magnetometry and polarised neutron 
diffraction. The updated phase diagram is substantially more complicated than previously re¬ 
ported and can be understood in terms of competing interactions with loops of spins, trimers and 
decagons, in addition to competition and interplay between antiferromagnetic, incommensnrate 
and ferromagnetic order. Several additional distinct phase bonndaries are presented. The phase 
diagram centers around a multiphase convergence to a single point at 0.9 T and ~ 0.35 K, below 
which, in temperature, a very narrow magnetically disordered region exists. These data illustrate 
the richness and diversity that arises from frnstrated exchange on the 3 dimensional hyperkagome 
lattice. 


1 



I. INTRODUCTION 


Frustration is ubiquitous in nature and drives the physical behaviour in materials ranging 
from liquid crystals and polymers to compounds with localised magnetic momentJ^. The 
study of magnetic frustration, in particular, is proving very fruitful in the development of a 
more general understanding of frustrated phenomena. Magnetic frustration is governed by 
the connectivity of the many degenerate spin conhgurations in its ground state manifold 
with geometric frustration borne out of the crystal structure topology. The multiplicity 
of ground states at the lowest temperatures, T —?• 0 K, can lead to persistent dynamic 
magnetic spin^, a most exotic state. 

The rare earth garnet Gd 3 Ga 50 i 2 (GGG) is unique since it offers a rare opportunity to 
study frustration on a double hyperkagome structure, a 3D kagome lattice of interconnected 
triangles. GGG is regarded as the archetypal frustrated compound since, unlike its many 
counterparts, it does not revert to an ordered state via an ’’order by disorder” transitiorP. 
Indeed there is no sign of long range order in GGG down to 0.025 K despite a Gurie-Weiss 
temperature of 9cw ~ -2 indicating strongly frustrated antiferromagnetic (AF) 

correlations. In contrast, the rare earth counterparts Dy 3 Ga 50 i 2 and Ho 3 Ga 50 i 2 order 
at relatively high temperatures thereby highlighting the unique nature of the magnetic 
exchange interactions in GGGP. 

The leading magnetic interactions in GGG are the near neighbor, Ji — -0.107 K®, and 
the dipole exchange interactions, D — 0.0457 Adjacent triangle and sublattice exchange 
interactions, J 2 and J 3 respectively, are an order of magnitude smaller, J 2 ~ -0.005 and J 3 
~ O.OIOK!^. In GGG the magnetic Gd^+ spins [S = 7 / 2 ) are considered as Heisenberg 
spins with single ion anisotropy of less than 0.04 K. However the non-negligible dipole 
exchange and the local crystal field environment could lead to anisotropjP^. The inherent 
spin Hamiltonian of the localised Gd^'*' moment on the hyperkagome lattice gives rise to 
magnetic short range order (SRO) below T r\j 3 KpHnUlIill Below about 0.018 K, a spin 
glass phase has been observed by ac susceptibility and specific heat at ambient field^ and 
verified by Petrenko et al^ using neutron scattering. Neutron scattering revealed sharper 
but not resolution limited magnetic diffraction peaks in addition to SRO for T < 0.14 hP. 


2 




These correlations are understood to result from the long range nature of the dipole exchange 
interactions perturbed by the very weak exchange interactions, J 2 and and may give 
rise to a quantum protectorate of ten ion spin clusterJ^. 

The first (H, T) phase diagram of GGG indicated a AF dome with strong anisotropy 
such that applying the field along the [1 0 0] creates an AF dome for T < 0.4 and 0.85 < 
PoH < 1.8 T while applying the field along [1 1 1] results in an AF dome for T < 0.3 K 
and 0.9 < /ioH <1.4 The anisotropy of the phase diagram was further probed via 
magnetic susceptibility, specific heat and single crystal neutron scatterin^^^^ and revealed 
that the phase diagram is much more complex than previously reported. The boundary 
between the disordered region above the spin glass phase and the long ranged AF dome 
showed considerable similarities with the melting curve of ^He indicating unusual magnetic 
behaviour in this regiorP^. In contrast, muon spin relaxation and Mossbauer spectroscopy 
indicate a temperature spin relaxation indicative of slow spin fluctuations down to 0.025 K 
and up to 1.8 ThSHZU. The contradiction between a long range ordered and dynamic magnetic 
state can be understood in terms of the time scales probed by the various techniques such 
that the system appears paramagnetic on the fast timescale of the muon but looks static 
within the energy resolution of the neutron. This has been observed in the frustrated magnet 
Tb 2 Ti 207 for which it was shown that there were many length and timescales to considei^^^. 

This study revisits the temperature and field dependent behaviour of GGG through 
macroscopic magnetic measurements of both powdered and single crystal samples and a 
polarised neutron scattering study of a powdered sample. In this study the field is applied 
along the [1 1 0] crystallographic direction to complement previous results. The phase 
diagram presents several additional phases for this direction of the applied magnetic field. 
These phases are also observed for the powder sample, but are broadened suggesting different 
characteristic field dependencies as a function of direction, as reported previously. Neutron 
powder diffraction correlate closely with the macroscopic measurements. 
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II. MAGNETIZATION MEASUREMENTS 


A. Experimental details 

Magnetization measurements were performed by the extraction method down to 0.08 
K and up to 8 T, using a superconducting quantum interference device magnetometer 
equipped with a miniature dilution refrigerator, developed at the Institut NeeP^. Two 
Gd 3 Ga 50 i 2 samples have been measured: i) a 8.52 mg single crystal cut from a commercial 
(10 X 10 X 0.4) mm^ substrate (Impex High Tech), measured along the [1 1 0] direction 
and ii) a 17.13 mg ^®°Gd isotope powdered sample. Neutron scattering experiments were 
performed on a powder sample of the same batch for a direct comparison. Magnetization 
measurements performed on the powder are in quantitative agreement with single crystal 
measurements. However, the features are broadened, certainly due to the distribution of 
crystallites and the presence of a small anisotropy 0. 04 KP). The results presented be¬ 
low are from the single crystal sample. The field was applied in the plane of the substrate, 
so demagnetization corrections are negligible. 

B. Phase diagram via magnetization measurements 

Low field magnetization measurements were performed by cooling down the sample from 
4.2 K to 0.08 K in an applied field of 10 mT. The susceptibility y (= M/H in low fields 
for which the magnetization is linear in field) as a function of temperature is shown in 
Fig. in The inset of Fig. [T] shows a monotonic decay of the inverse susceptibility with 
decreasing temperature. A linear behaviour is observed down to 1.5 K and assumes a 
Curie-Weiss law between 1.5 and 4.2 K that provides a Gurie constant of 8.13 emu/mol.K, 
consistent with the Gd^"*" spin moment, S = 7/2. The obtained Gurie-Weiss temperature 
9 is -1.97 K in agreement with previous resultd^^. Below 1 K, the plot deviates from the 
linear behaviour suggesting the development of correlations between the spins. In Figj^ 
isothermal magnetization measurements shows saturation around Mgat = 7 /is/Gd. Above 
2 K, no anomaly is seen in the magnetization curves. 

A more complex behavior occurs below 1 K. To analyse the M{H) curves as a function of 
temperature, dM/dH is plotted as a function of the applied field, see Fig. The obtained 
features show much correspondence with the observations reported by Schiffer et al^ in 
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FIG. 1: X vs r at 0.01 T with H || to [1 1 0] for 0.08 <T < 4.2 K. Inset: x ^ vs. T. 



iioH [T] 

FIG. 2: M vs for 0 < /ioH < 7.5 T with H || to [1 1 0] for several temperatures between 0.08 

and 4.2 K. Inset: zoom on the anomaly between 0 and 1.5 T at low temperature. 

heat capacity and susceptibility measurements. Fig. shows the main features with a 
broad maximum for H' ^ 1 T that develops below 1 K. It has previously been assigned 
to the quenching of the AF short range order (SRO). There is no significant change in the 
position or the width of this broad feature with decreasing temperature. Three well-defined 
peaks labeled Hi, H 2 and develop as the temperature is reduced below ~ 0.34 K in 
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addition to the broad feature centered around 1 T. The peaks at the lowest field [Hi = (IqH 
~ 0.65 T) and at the highest field (H 3 = [IqH ~ 1.25 T) correspond to the boundaries of the 
previously observed field induced AF long range ordei^. The temperature dependence of 
the peak positions are in good agreement with previous measurements, but their positions 
are slightly different possibly the result of anisotropy. H 2 is a very clear and previously 
unreported peak at around H 2 = HqH ~ 0.9 T. A weak feature was also observed in this field 
range by Schiffer et al^, but it was much more rounded and of smaller amplitude.These three 
peaks, Hi, H 2 and H 3 , grow on top of the broad SRO feature at H' but do not suppress it, 
confirming the coexistence between long-range and short-range ordering. It is worth noting 
that the newly reported H 2 peak is very sharp, similar to the Hi and H 3 peaks which were 
identified as phase boundaries of the field induced ordered phas#*^ thereby indicating that 
H 2 is the signature of a phase boundary. With increasing temperature the positions of the 
three peaks and of the broad feature converge towards a single point around T* ~ 0.35 K 
at the magnetic field fioH* ~ 0.9 T. To get a deeper understanding of this feature, we 
performed magnetization measurements as a function of temperature at a constant field. 



m 


FIG. 3: dM/nodH vs. hqH for a range of applied fields 0 < hqH < 8 T with H || to [110]. The 
curves for different temperatures are separated by 1 ps for clarity. The intensity has been fitted 
with a Gaussian function for the broad feature and three Lorentzian functions for the three peaks. 

At low fields, below 50 mT, the magnetization shows a freezing below 200 mK in the 
zero field cooled - field cooled (ZFC-FC) measurements, as previously reportecP^. This 
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FIG. 4: X vs. T measured with a ZFC-FC procedure for a rauge of applied fields 0.4 < /noH < 1.3 T 
aloug the [1 1 0] directiou. 


freeziug was also observed iu our ac susceptibility measuremeuts. Above 50 mK, this ZFC- 
FC irreversibility disappears. Iu the followiug we focus ou larger applied fields to further 
explore the phase diagram iu comparisou with our M versus H measuremeuts. Fig. 
shows ZFC-FC temperature depeudeut maguetizatiou measuremeuts for a rauge of applied 
maguetic fields, Fig. ia): 0.4 < iJoH < 0.7 T , Fig. IJb): 0.8 < fioH < 0.9 T and 
Fig. 1^ (c):1.0 < fioH < 1.3 T. The magnetic field dependence can be subdivided into 
five distinct behaviours. First: (IqH < 0.4 T. ZFC and FC magnetizations overlay and 
decay monotonically with increasing temperature. No magnetic ordering is observed in 
this field range. Second: For 0.5 < hqH < 0.7 T, two distinct features are observed: a 
broad maximum around Ti — 0.6 K and a change in slope around T 2 = 0.2 K. Ti shifts 
to lower temperatures when the magnetic field is increased. Third: For 0.8 < fioH < 0.9 
T, the two features overlap so that no clear maximum emerges. A splitting in the ZFC- 
FC measurements occurs at fioH* = 0.9 T below T ~ 0.2 K. This ZFC-FC hysteresis is 
quantitatively reproducible and is present in a very narrow field region around 0.9 T. The 
ZFC-FC splitting concerns 7% of the magnetization. Fourth: For 0.9 < hqH < 1.2 T, the 
curves recovers the shape of the 0.5 < (IqH < 0.7 T region. Fifth: For [IqH > 1.3 T, the 
magnetization decreases continuously with increasing temperature. 

An updated phase diagram can be constructed from these measurements and is shown in 
Fig. [^and Table IIB, High temperatures, Ti (orange circles) above H' (violet squares) can 
be associated with the SRO features previously observed in specific heat and susceptibility 
measurements by Schiffer et alW^. However, the low temperature change in slope around T 2 
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FIG. 5: The updated phase diagram of Gd 3 Ga 50 i 2 determined from macroscopic magnetization 
measurements and neutron powder diffraction. The continuous white lines show regions measured 
in this study using neutron diffraction. 

(red squares) indicates an evolution of the SRO state, which thus can be separated into two 
regions (A-E) separated at the onset of T 2 . Hi and H 3 (blue triangles) match the AF field 
induced region previously observed. The new observed peak in M vs fioH measurements 
at H 2 = fJ^oH = 0.9 T (violet triangles) divides the AF field induced region into two phases 
(B-C). At this (B-C) phase boundary, [IqB. ~ 0.9 T, T < 0.2 K, a small region is reported 
in which a ZFC-FC discrepancy is observed. Interestingly this phase boundary corresponds 
to the development of field induced phases and the collapse of SRO more explicitly probed 
using neutron scattering. As such it makes /iqH = 0.9 T, T ~ 0.35 K a specific point in 
the phase diagram of GGG. In the following section these results are compared with the 
coexistence of short range, incommensurate and AF phases as observed by Petrenko et. 
and extended in this work using neutron diffraction. 




1 1 

A dM/|JodH(H) (LRO - upper peak) 
A dM/ModH(H) (LRO - middle peak) 
A dM/podH(H) (LRO - lower peak) , 

• dM/podH(H) (SRO - H’) 

♦ M/|JoH(T) (T1) 

■ M/^JoH(T) (onset of T2) 

I ^ Neutron diffraction ; 

0 

ODo „ 

F 

0. 





c 

AAi 


i • ZhU-h 

Ref. n 

■0 discrepanc 

71 


P.« ... 







0 ^ 


0 

1 

B 


► 

w ^ 

0 

0 


A 


E 

. 0 . 

0 

0 

0 

0 


1 

i_i 



i J 

























Phase 

T [K] fioHlT] 

M Signature 

NS Signature 

A: (T < 0.25, noH < 0.65]) Slope change (T 2 ) 

• Broad feature at 0.8 

• Broad feature at 0.8 

• Development of IC and AF peaks 

E; ([0.25; 1 ], fioH < 0.9 ) 

Broad maximum (Ti) 

A Loss of AF (1.14 A-i) & IC (1.08 A"!) 

A Decrease of SRO 

B: ( T < 0.35, [0.65; 0.9 ] ) Three peaks in 

• IC peak ^[2 10] 

• Loss of IC peaks between [1 1 0] ^[| 1 |] 

C: ( T < 0.35, [0.65; 0.9 ]) 

Irreversible ZFC-FC 

A IC peak between [2 0 0] and [2 1 0] fixed in Q 

A Scattering at Q = 1.69 A“^ 

A SRO more correlated 

D: ( T < 0.2, /roH > 1.2 ) 

No signature 

• Development of scattering at 0.8 A“^ 

• Loss of first SRO peak 

• Appearance of IC peak at Q = 0.80(2) A“^ 

• Loss of Bragg peak at 1.69 A“^ 

F: ( [0.2;1], > 0.9 ) 

No signature 

A Loss of [2 1 0] scattering 

A Weak variation of [2 0 0] Bragg peak intensity 


TABLE I: Outline of phase diagram presented in Fig. SRO = short range order, M = magne¬ 
tization, NS = neutron scattering. 

III. NEUTRON DIFFRACTION 

This study focusses on the regions of the phase diagram encompassing 0.175 < T < 3 K 
(HoH = 0 T) and 0.06 < T < 3 K ( 0 < hqH < 2 T). The nnnsual scattering observed for 
T < 0.14 K in zero field with longer ranged order snperposed on short range correlations is 
not considered. Neutron scattering corroborates the complicated phase diagram shown in 
Fig. [^with strong consistency between magnetization and neutron diffraction measurements. 

Petrenko et al^ conclnded that competing magnetic interactions result in competition 
between various ground states and prevents GGG from ordering in zero magnetic field while 
applying a small magnetic field crystallises the magnetic state into ordered components with 
a range of ferromagnetic (FM), AF and incommensurate (IC) propagation vectors. 
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The present work provides a greater insight into the interplay between different magnetic 
orders. Particular emphasis is placed on the perturbation of the short range order with 
respect to the other magnetically ordered states. Short range magnetic order pervades 
the {H, T) phase in which uniquely long range magnetic correlations have previously been 
assigned. The interplay between short range, FM, AF and IC orders are correlated with the 
phase diagram outlined in Fig. [^and with the main results condensed in Table IIB These 
datasets highlight why a refinement of the magnetic structure with a single wavevector has 
eluded previous studies. 


A. Experimental details 

Neutron scattering experiments were performed on the D7 diffuse scattering spectrometer 
and the cold chopper spectrometer IN5 at the Institut Laue Langevin (ILL), Grenoble. The 
sample has previous been used in the work of Petrenko et al^ and Been et al^ containing 
99.98% of the non-absorbing isotope ^®°Gd. The sample was covered with Iso-Propanol (or 
2-Propanol) 99% Deuterium that freezes the crystallites into place without any substantial 
contribution to the scattering. 

The zero field data, measured on IN5, shows the scattering within the elastic resolution 
of the instrument, 80 /reV full width at half maximum (FWHM). The energy resolution is 
determined using a standard incoherent scatterer, see section [Til C 1[ The INS data set could 
not be integrated in energy due to spurious scattering in parts of the inelastic spectrum. 

The field dependent data has been measured using D7 with E/ = 3.55 meV. The wide 
angular range available on D7 enables simultaneous determination of short and long ranged 
scattering and elucidate their interplay. Field dependent measurements up to 2.5 T were 
performed at 0.05 K with a temperature dependence measured, up to 0.8 K, at 0.8 and 
1.2 T, as depicted by the white continuous lines on Figj^ Uniaxial polarization analysis 
performed under applied field on D7 provides two scattering cross sections, spin-flip and 
non-spin flip scattering. The spin-flip scattering is entirely magnetic in origin since the spin 
incoherent cross-section for GGG is negligible. The non spin-flip scattering contains nuclear 
and magnetic components. A separation of the spin-flip from the non-spin flip contributions 
allows the nuclear scattering to be isolated while the spin-flip contributions provide the 
magnetic cross-section^. The data have been corrected for detector and polarization an- 
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alyzer efficiencies using standard samples of vanadium and amorphous silica respectiveljPSl. 
All data are corrected for background contributions by subtracting the scattering from an 
empty sample can under equivalent conditions. The D7 data set is integrated in energy. 
Nevertheless the data sets from D7 and IN5 can be compared since previous work show that 
the majority of the scattering is elastic (82%) and the features of interest for this study he 
solely in the elastic This is discussed further in the appendix. 


B. Ambient field behaviour 


Region A, depicted in the phase diagram of Fig. has generally been considered to 
consist of short range magnetic correlations. Petrenko et determined, using Monte 
Carlo (MC) methods, near neighbor interactions, Ji = 0.107 K, as the dominant magnetic 
exchange term. Upon closer inspection of the scattering profiles, see Fig. |^a), it can 
be seen that a model with only near-neighbor exchange cannot give rise to the scattering 
profile observed. Instead, it is proposed that a self organised super-lattice unit, a ten ion 
spin cluster, is required to account for some of the features, see Fig. |^b). The arguments 
behind such a model are presented in the next few paragraphs. 

Short range magnetic correlations in a powder sample scatter according to the expression 



(So • S„) 

S{S + 1) 


N, 


sm{QR„ 

QRn 


-F{Q), 


( 1 ) 


where Sq and Sn are the spin magnitudes of the central atom and the shell atom, R^ and 
Nn are the radii and coordination numbers of the nearest neighbor shell respective!}^ 
and F{Q) is the magnetic formfactoi®3SI j^a) shows the expected magnetic scattering 
profile for short range order, using eqn. [!} in comparison to the magnetic scattering profile 
measured at T = 0.25 K (0 T). Interestingly, the short range order profile does not reproduce 
all the features expected, in particular, the broad feature centered on Q = 1.8 A“^. In order 
to reproduce scattering at Q = 1.8 A magnetic exchange between ten ion spin clusters is 
considered. 

In recent years there has been much excitement surrounding emergent behaviour in geo¬ 
metrically complex materials. The kagome lattice has provided such an emergent structure 
in the form of correlated hexagon loops of spin^. It has been suggested that the first and 
second order zero energy modes on the hyperkagome structure are dynamic loops of spins 
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that include triangular loops of spins, trimers, and ten ion spin clusters, decagond^^®. The 
interpretation from recent optical hole-burning experiments supports the existence of ten ion 
spin clusters in the low temperature spin glass phase of GGGI^. Fig. [^b) shows a snapshot, 
in time, of possible spin directions involved in such a ten ion spin cluster for the [1 0 0]/[0 1 0] 
crystallographic directions. The loops are decoupled from their environments and suggested 
to be quantum protectorate^^®. The magnetic exchange for trimers can be modelled using 
eqn. [TJwith only near neighbor exchange and the requirement that Heisenberg spins on a tri¬ 
angular lattice '^triangle S = 0 resulting in near neighbour spins 120" relative to each other. 
The magnetic neutron scattering profile for a trimer, see Fig. [^a), shows a maximum in 
scattering intensity at the same Q position as near neighbor magnetic correlations, Q ~ 1.1 
A“^. The spins on a decagon loop also obey the Heisenberg near neighbor spin requirement 
yet the spins are correlated throughout the ten ion spin cluster. The exchange between ten 
spin ion clusters can be described using eqn. [T] to consider short range interactions between 
adjacent clusters. The coordination number of a ten ion spin cluster is 14 with an average 
unit cell size of 13.12 A. The magnetic scattering from a ten ion spin clusters, see Fig. |^b), 
provides many extra features in common with the powder diffraction data of GGG at 0.25 K 
and 0 T. 

The models show that using a combination of near neighbour and exchange between ten 
spin ion clusters can explain some experimental features. First, around Q ~ 0.3 A“^ there 
is a weak feature consistent with the position of the first peak of the ten ion spin cluster, 
see inset Fig. [^a). Second, the main peak of the near neighbor model is found at the same 
position as the main peak in the data, however this peak is strongly asymmetric and this 
can be ascribed to scattering from ten ion spin clusters which provides two features either 
side of the main peak position, see dashed lines in Fig. |^a). Thirdly, a strong scattering 
feature at 1.8 A“^ in the ten ion spin cluster model is consistent with the data. Interestingly, 
the magnetic correlation length ^ = 2 tt/AQ of the SRO peak at 1.8 A“^ corresponds to an 
average unit cell size of the ten ion spin clusters, ^ = 13.8±2 A. A combination of these two 
models captures the position of the scattering if not perfectly the respective weights of the 
scattering. Nevertheless it provides a strong indication of unusual correlations such as the 
coexistence of near neighbor exchange with a ten ion spin cluster. This dataset therefore 
captures the first two soft modes of spins on hyperkagome structure, dynamic trimer and 
decagon spin loops as a feature of region A. 


12 





(b) 


FIG. 6: (a) Neutron diffraction of Gd 3 Ga 50 i 2 at 0.25 K, 0 T with a near neighbor (NN) spin 
model (red line -) and exchange between ten ion spin clusters (black -.-) for comparison. Inset: 
Zoom of the magnetic scattering around Q = 0.3 A~^. (b) Snapshot of ions involved in the ten ion 
spin cluster for the [1 0 0]/[0 1 0] crystallographic direction. The blue rectangle is the structural 
unit cell. The light green and light blue loops show two domains. The arrows denote the relative 
spin directions of the ions. 

C. Magnetic behaviour under applied magnetic fields. 

1. Field dependence (T = 0.05 K) 

An overview of the field dependence of the magnetic scattering is shown in Fig. 
and reveals great complexity. Peaks corresponding to AF peak positions are marked with 
continuous blue lines while FM peak positions are marked by dashed black lines. The 
phase boundaries reported on Fig. [^are reproduced on Fig. Variations in the magnetic 
scattering profiles can be accurately mapped onto the phase diagram as described in Table 
ng A closer look at the individual scattering profiles for each applied magnetic field 
provides more detailed information. 


Fig. i a) shows magnetic scattering profiles in region A —)• H for 0, 0.2 and 0.5 T 
(T < 0.25 K). The scattering profile at 0.25 K, 0 T is the elastic scattering within the 



13 


















0.6 0.8 1 1.2 1.4 1.6 1.8 

Wavevector Transfer [ A ~ ^ ] 


FIG. 7: Overview of the field dependence of the magnetic neutron scattering at 0.05 K. Peaks 
corresponding to AF peak positions are marked with full vertical blue lines while FM peak positions 
are marked by dashed vertical black lines. The phase boundaries found in Fig. [^are presented by 
dashed horizontal lines. 

energy resolution, 80 /ieV of INS, further information is provided in the appendix.Upon the 
application of only 0.2 T the short ranged ordered peaks are drastically reduced. This is 
particularly true for the first SRO peak at Q ^ 1.1 A“^. Interestingly, the second SRO 
peak, assigned to scattering from a decagon loop in the previous section, is more robust 
and remains unaffected up to 0.5 T. This would indicate that decagons are more resilient to 
perturbation than the trimer counterparts upon the application of a magnetic field. Longer 
ranged correlations are developed by 0.2 T and are resolution limited at 0.5 T, see Fig. 

Fig. a a) shows the complete field dependence of the magnetic diffraction patterns 
at base temperature, 0.05 K. FM correlations are not considered in this work since it 
is impossible to disentangle the effect of local magnetic exchange interactions and the 
increase in scattering at the nuclear Bragg peak positions due to the application of a 
magnetic field on a powder. The data are offset for clarity. The continuous vertical lines 
correspond to positions in reciprocal space where AF order is expected while dashed 
vertical lines are positions corresponding to nuclear/FM order. In region A, 0 < /ioH < 
0.65 T, the IC features observed between the [1 1 0] and [| 1 f] Bragg peak positions 
and for Q < 1.2 are very fluid with respect to their absolute position in reciprocal 

space. This is particularly true for the broad features around 0.6 - 0.8 A~^ at 0.2 T 
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FIG. 8: (a)Field dependence of magnetic neutron diffraction at 0.25 K & 0 T, 0.05 K & 0.2 T and 
0.05 K & 0.5 T. 


T = 0.05 K 



(b) 


G 
G 

■£ 
^10- 


Magnetic 


[3 1 1] 
Forbidden 

i ic 

o 

oo. 


^ goooooo 


o 


o 


(D J 


Nuclear 


q^OoOoOQQqoOOOOO 


OOqOO 

oo^oop 


1.6 1.65 1.7 1.8 

Wavevector Transfer [A'*] 


1.9 


FIG. 9: (a) Field dependence of the magnetic scattering at 0.05 K. The continuous blue vertical 
lines correspond to order that is commensurate with the crystal structure and could represent AF 
order. Dashed black vertical lines are nuclear/ FM positions. The vertical arrows denote scattering 
at the forbidden nuclear scattering position [3 1 1] and an IC position. The data are offset for 
clarity, (b) Scattering profiles of the magnetic and nuclear contributions at 0.05 K and 1 T zoomed 
in on the region of the scattering at 1.69 and 1.72 A“^. The nuclear profile shows no equivalent 
scattering. 
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that have shifted to between 0.7 - 0.9 A~^ by 0.5 T and are not observed beyond. The 
boundary between phases A and B fixes the Bragg positions of various AF and IC Bragg 
peaks. An example is the [2 1 0] Bragg peak position which appears to move from an 
IC position in phase A to the [2 1 0] Bragg peak position in phase B. Additionally, the 
IC scattering between [2 0 0] and [2 1 0] is also fixed into position as the A-B phase 
boundary in crossed. Upon the establishment of phase C, 0.9 < poH < 1.2 T, magnetic 
scattering for Q < 0.8 is diminished and magnetic scattering appears at Q = 1.69(2) 
A~^. Interestingly this scattering corresponds to the [3 1 1] Bragg peak position, forbidden 
in neutron scattering for this crystallographic symmetry. This scattering could then be 
the result of either a magnetic phase change or a signature of magneto elastic coupling 
(MEC). MEC can be a very weak effect, possibly the result of a slight rotation of oxygen 
around the dodecahedral Gd^+ site. However, MEC should also give rise to scattering 
in the nuclear channel however neutron scattering might not be sufficiently sensitive to 
reveal the nuclear scattering part. Indeed there is no nuclear scattering peak at the [3 
1 1] Bragg peak position when comparing the magnetic scattering cross section to the 
nuclear scattering cross section, see Fig. [^b), so the origin of this scattering remains unclear. 

Region D is consistent with the complete loss of the first SRO peak, the appearance of 
a IC Bragg peak at Q = 0.80(2) and the loss of the Bragg peak at Q = 1.69(2) A“^. 
Although the first SRO peak has disappeared in this region, the second SRO peak, at Q = 
1.8 A“\ remains strong up to and beyond 2.0 T. Above 1.5 T, all AF peaks are suppressed 
except the AF peak at 1.01 A ([2 0 0]) which persists. 

An overview of the field dependent behaviour of the integrated intensities of some of the 
magnetic scattering is given in Fig. It is particularly interesting to note the interplay 
between the various magnetic contributions such that the AF scattering at Q = 1.14 A“^, 
([2 1 0]), is most intense when the AF scattering at Q = 1.01 A“^ ([2 0 0]) is much reduced 
only for the Q = 1.14 A“^ peak to restrengthen when the Q = 1.01 A“^ peak diminishes, 
region D. In addition, the field dependence of the IC peak at 1.08 A“^ follows the integrated 
intensities of the AF peak at Q = 1.14 A“^. 
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FIG. 10: Field dependent integrated intensities of AF, IC and short range order. The data is 
normalised for ease of viewing and the lines through the data are to guide the eye. 

2. Temperature dependence at 0.8 T (T = 0.05 K) 


The magnetization data implies a convergence phases to a unique point around 0.9 T and 
0.35 K. Temperature dependent powder diffraction has been measured at 0.8 T to determine 
the unique nature of this region, see Fig. [IT|[a,b), with particular emphasis on the region 
0.2 < T < 0.4 K. These data show a very distinct transition at 0.3 K with both AF and 
IC orders disappearing beyond 0.3 K, see Fig. inic). The line separating regions A —>■ E 
is observed in the neutron scattering profiles. The magnetic scattering correlation of the 
first SRO peak broadens beyond 0.2 K from 13.1(3) A at the lowest temperatures to 9.2(2) 


A at 0.8 K, see inset of Fig. |ll[c). It is difficult to extract ^ for the second short range 
ordered peak due to its weaker nature and limited Q-range. The magnetic scattering at the 
IC position, Q = 1.73 A“^, and the forbidden nuclear Bragg peak position Q = 1.69 A“^ 
disappear between 0.2 and 0.3 K, see the black arrows on Fig. inib)- These data do indeed 
indicate that there is a convergence of phases around 0.9 T and 0.35 K as indicated in Fig. 
however more detailed studies are needed to elucidate the precise details of this region. 


3. Temperature dependence at 1.2 T 

The neutron scattering data of the 1.2 T temperature dependence. Fig. [T^ traverses 
regions C ^ D ^ F. The IC Peak at 1.08 A disappears between 0.2<r<0.3K, black 
arrow in Fig. [T^b), region C ^ D. The peak at the position of the forbidden nuclear Bragg 
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FIG. 11: (a) Temperature dependence of the magnetic scattering at 0.8 T. The data are offset 

for clarity. The blue dashed vertical lines corresponding to AF order. Dashed vertical lines are 
nuclear/ F positions, (b) Focus on the region 0.2< T < 0.4 K (0.8 T). The black arrows emphasise 
the greatest variations in this temperature range, (c) Temperature dependence of the integrated 
intensities of the AF (Q = 1.14 ) and IC (Q = 1.08 ) clearly showing the phase transition 

A ^ E defined in Fig. Inset: Temperature dependence of ^ for the first SRO peak. 

peak, Q = 1.69 A“^, in addition to the AF Bragg peak at Q = 1.14 A ^ disappear between 
0.3 < T < 0.4 K, region D ^ F. Interestingly the peak at Q = 1.01 A“^, a peak that can be 
attributed to AF order, is only slightly perturbed in the 0.3 < T < 0.4 K region. Fig. [T^c) 
shows the variation in scattering prohles across the F —>■ E transition. A slight variation in 
intensity is observed but no clear transition can be observed via neutron diffraction. Single 
crystal measurements are required to elucidate the details of the variation in scattering to 
determine whether the regions E and F are distinct phases. 


IV. CONCLUSION 

An investigation of the applied magnetic field and temperature dependent behaviour of 
the archetypal magnetically frustrated compound Gd 3 Ga 50 i 2 is presented. Single crystal 
and powder magnetic susceptibility in conjunction with polarised neutron diffraction reveal 
a unified picture of the phase diagram with close correlation between magnetic susceptibility 
measurements with the field applied along the [1 1 0] crystalline direction and neutron 
scattering profiles of a powdered sample. 
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FIG. 12: (a)Temperature dependence of Gd 3 Ga 50 i 2 at 1.2 T. The data are offset for clarity. The 
blue dashed vertical lines corresponding to AF order. Black dashed vertical lines are nuclear/ 
FM positions, (b) Integrated intensity of AF and IC Bragg peaks. Lines are a guide to the eye. 
The arrows show regions are greatest variation in this temperature range, (c) Neutron scattering 
profiles across the phase boundary E —> F. No clear transition is apparent. 

Several extra phase boundaries are required to correctly describe the (H,T) phase diagram 
of GGG with H II [1 1 0] . There is strong evidence that trimerised and decagon loops of 
spins coexist at low temperatures and low magnetic fields. These emergent loops of spins are 
strongly affected by the long range ordered components yet coexist with incommensurate 
and AF order up to 1.3 T for spin trimers and beyond 2.0 T for decagon spin structures. 
Magnetization measurements allude to a multiphase convergence around 0.9 T and 0.35 K. 
Interestingly, there is a strong ZFG-FG discrepancy at 0.9 T and T < 0.2 K that is very 
well defined in applied field indicative of a further glassy phase. The close link between 
susceptibility and neutron powder scattering thus indicates that the behaviours observed 
exist throughout the crystal, originates from the spin Hamiltonian on the hyperkagome 
lattice and are an intrinsic feature of frustration on hyperkagome structures. 
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